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Genetically Diverse Coronaviruses in Captive Bird 
Populations in a Brazilian Zoological Park 


Abstract 

This study aimed to investigate the occurrence of coronaviruses (CoVs) in captive birds placed inside a zoo¬ 
logical park in Brazil. The role of captive birds in the epidemiology of Co Vs in the tropics is poorly understood. 
A total of 25 (ft = 25) different species were tested for viral RNA using individual fecal samples collected from 
healthy birds. Reverse transcription-polymerase chain reaction targeting the 3' untranslated region was used to 
detect CoV RNA, and positive samples were submitted for sequence analysis. The phylogenetic search revealed 
nine mutations in the black shouldered peafowl (.Pavus cristatus) CoV sequence, which clustered separately from 
samples previously described in England. This is the first report on the detection of the CoV genome in captive 
birds in Brazil. 
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Introduction 

T he outbreaks of highly pathogenic avian influenza 
viruses around the globe in the past few years have been 
devastating for the countries concerned and have created 
fear in the general population. Similarly, in 2002 the world 
witnessed the emergence of severe acute respiratory syn¬ 
drome in Guangdong Province, China, caused by a previ¬ 
ously unknown coronavirus (CoV) (Liu 2003). It is likely that 
additional, yet to be identified CoVs are circulating in do¬ 
mestic birds, especially in tropical countries. The three ex¬ 
isting groups of Co Vs were initially distinguished based on 
the lack of an antigenic relationship between them (Masters 
2006). Group 3 CoVs, which are genetically similar to in¬ 
fectious bronchitis virus (IBV), have been detected in healthy 
galliform and nongalliform birds (Cavanagh et al. 2002). 
Recently, populations of wild birds in northern England 
were found to have IBV-like RNA in their feces (Hughes 
et al. 2009). Currently, little information is available about 
the role of South American bird populations as reservoirs or 
long-distance vectors of CoV-like viruses. 

Viral diseases can cause substantial mortality in both wild 
and captive bird populations. Brazil is a country with many 
migrating wild birds, which is likely to facilitate disease 
transmission from migratory to nonmigratory flocks in zoo¬ 
logical parks or on farms. As viral detection and identification 
techniques become increasingly sophisticated, and as the 


study of wildlife diseases increases, it is likely that new 
viruses will be discovered and new diseases will be encoun¬ 
tered (Hughes et al. 2009). Therefore, studies of the viral dis¬ 
eases of wild and captive birds are necessary. The latest 
techniques of molecular epidemiology can help scientists 
understand how viruses move from species to species and 
from region to region (Hughes et al. 2009, Jackwood et al. 
2010 ). 

The present study demonstrates the molecular detection 
and characterization of viruses from captive exotic bird pop¬ 
ulations sampled in Sao Paulo State Zoological Park in Brazil, 
which appear to be related to group 3 CoVs. 

Materials and Methods 

Study design 

Serial cross-sectional surveys of captive wild bird pop¬ 
ulations, maintained in a zoological park in the city of 
Bauru (N22°10 / S49°00 / W), were undertaken in July 2009. 
Feces were collected from 25 different species of healthy 
birds (Table 1), placed in virus transport media (Eagle's 
minimum essential medium supplemented with 10% fetal 
calf serum, penicillin, and streptomycin suspension, plus 
amphotericin B), and frozen at —80°C for individually 
analysis later. This work was authorized by the Animal 
Ethics Committee of the Sao Paulo State University (pro¬ 
tocol number 2008/09878-09). 


Laboratory of Animal Virology, UNESP, School of Veterinary Medicine, University of Sao Paulo State, Aragatuba, Brazil. 
2 Coronavirus Research Group, School of Veterinary Medicine, University of Sao Paulo, Sao Paulo, Brazil. 
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Table 1. Captive Bird Species from Whom Fecal Samples Were Collected and Screened for Coronavirus RNA 


Common name 

Latin name 

Age 

Positive/total samples 

Chestnut-bellied heron 

Agamia agami 

Ad. a 

1/20 

Egyptian goose 

Aiopochen aegyptiacus 

Ad. 

0/20 

Amazon parrot 

Amazona estiva 

Juv b 

0/5 

Blue-front parrot 

Amazona vinacea 

Ad. 

0/5 

Bar-headed goose 

Anser indicus 

Juv. 

0/5 

Emperor penguin 

Aptenodytes patagonicus 

Ad. 

1/2 

Cattle egret 

Bubulcus ibis 

Juv. 

0/0 

Red-legged seriema 

Cariama cristata 

Ad. 

0/4 

Kelp gull 

Larus atianticus 

Ad. 

0/2 

Silver pheasant 

Lophura nycthemera 

Ad. 

0/1 

Tawny-bellied screech-owl 

Megascops watsonii 

Ad. 

0/12 

Red crested cardinal 

Paroaria coronata 

Juv. 

0/2 

Black shouldered peafowl 

Pavus cristatus 

Ad. 

2/5 

Indo-Chinese Green peafowl 

Pavus muticus 

Ad. 

0/14 

Ring necked pheasant 

Pavus torquatus 

Ad. 

0/1 

Dusky legged guan 

Penelope obscura 

Ad. 

0/1 

Rusted margined guan 

Penelope supercilliaris 

Ad. 

1/2 

Pheasant 

Phasianus colchicus 

Ad. 

0/12 

Greater flamingo 

Phoenicopterus rubber 

Ad. 

0/9 

Black-fronted pippin-guan 

Pipile jacutinga 

Juv. 

0/8 

Caraca cheryway 

Plyborus plancus 

Ad. 

0/23 

American purple gallinule 

Porhyrula matinica 

Ad. 

2/3 

Saffron toucanet 

Pteroglossus bailloni 

Ad. 

0/16 

Lettered aracari 

Pteropglossus inscriptus 

Juv. 

0/12 

Toco toucan 

Ramphastos toco 

Ad. 

0/10 

Channel-billed toucan 

Ramphastos vitelinus 

Juv. 

0/10 

Greater rhea 

Rhea americana 

Ad. 

0/10 

Roadside hawks 

Rupornis magnirostris 

Ad. 

1/10 


a Ad., adult (>2 years). 
b Juv., juvenile (<2 years). 


Viral RNA extraction 

Two grams of feces was homogenized in two volumes of 
Eagle's essential medium (Invitrogen™; pH 7.4) and clarified 
by centrifugation at 3000 g for 20 min. The supernatant was 
filtered once through a 0.75-pm paper filter (Millipore™) and 
then twice through a 0.25-^m syringe filter (Corning™) as 
described previously (Teixeira et al. 2007). These suspensions 
were heated at 100°C for 30 min in a water bath before total 
RNA extraction was performed. Total RNA was extracted 
using Trizol® standard protocol, based on guanidinium iso¬ 
thiocyanate and acid-phenol with some modifications (Teix¬ 
eira et al. 2007). We mixed 200 /iL of clinical suspension with 
500 pL of Trizol reagent and incubated for 10 min at room 
temperature. After addition of 200 /iL of chloroform, the re¬ 
action was mixed vigorously for 10 s and centrifuged at 13,000 g 
for 10 min. The upper aqueous phase was mixed with an equal 
volume of cold isopropanol and incubated on ice for 10 min. 
The total RNA precipitate was pelleted by centrifugation at 
13,000 g for 10 min and washed with ethanol. The RNA was 
dissolved in 30 pL of diethylpyrocarbonate-treated, sterile, 
double-distilled water and stored at —20°C. 

Reverse transcription-polymerase chain reaction 

The primer combination untranslated region (UTR) 11 — / 
UTR41+, which corresponds to the conserved region of the 
IBV 3' UTR and produces a very sensitive reverse transcription- 
polymerase chain reaction (RT-PCR) result (Adzar et al. 1997, 
Culver et al. 2006), was used. Before RT-PCR was performed. 


RNA was heated for 3 min at 100°C, followed by 2 min at 
72°C, and finally for 2 min at 50°C. Optimization of a one-step 
RT-PCR requires compromises in the concentrations of PCR 
reagents, particularly the primers used. The RT-PCR was per¬ 
formed following the instructions in the One-Step Superscript 
HI® pqvpcp commercial kit. Briefly, 2 pL of the first-strand 
product was used as a template for amplification in a 50 /iL 
reaction mixture containing 20pmol of primers and 1.6 mM 
MgCl 2 . All PCR conditions were the same as described above, 
but the amount of RT-Taq DNA polymerase platinum enzymes 
was calculated to 2.5 units/50 pL reaction mixture. A total of 
10 of PCR products was electrophoresed at 100 V for lh 
using a 1.5% agarose gel in Tris-borate EDTA buffer and vi¬ 
sualized with ethidium bromide staining and an ultraviolet 
transluminator. Gel images were captured using a Kodak 
DC290 digital camera and ADOBE 6.0 software. 

DNA sequencing and phylogenetic analysis 

The PCR products were purified using a commercial pu¬ 
rification kit (QIAquick PCR purification kit; Qiagen) ac¬ 
cording to the manufacturer's instructions and sequenced 
directly with an ABI 3730x1 DNA analyzer (Applied Biosys¬ 
tems). Nucleotide sequences derived from this study have 
been deposited in the GenBank sequence database under 
the following accession numbers: GU323504, Agamia/BR/ 
UNESP/2009; GU323505, APatagonius/BR/UNESP/2009; 
GU323506, Pavus/BR/UNESP/2009; GU323507, Pavus/2/ 
BR/UNESP/2009; GU323508, Penelope/BR/UNESP/2009; 
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GU323504 (Agamia/BR/UN ESP/2GQ9)* 

GU 3235 10 (M a rti n i c a /B R /U N ES P/2009) * 

__ GU323509 (Martinica/BR/UNESP/2009)* 

GU32351 1 (Ru porn is/BR/UN ES P/2009)* 


-GU323508 (Penelope/BR/UNESP/2009)* 

AJ27S335 (IBVD207) 


so 


GU323506 (Pavus/BR/UNESP/20G9) + 
GU323505 (A P ata gon i u $/ B R/U N E S P/2009) * 

AJ 27 B 33 S (IBVHVI-140) 

AJ27S337 (IBVHV10) 


AJ3i 0642 (T 0 oV turk ey AJ K/412/00)* 

-‘Ajelaseo (Ph/UK/602/95) 


69 


AJ 619579 (Rh AJ K/438/94) 

= GU323507 (Pavus/2/BR/UNESP/20Q9)* 


]--—i 

101 

FIG. 1. Minimum-evolution tree of corona viruses (CoVs) based on a 146-nt fragment of the 3' untranslated region from 
group 3 CoVs. Evolutionary distances were computed using the Tamura-Nei method and are expressed as the number of 
base substitutions per site. CoVs detected in captive wild birds by this study are denoted with an asterisk. Previously 
published CoV sequences from different sources were included for comparative purposes. GenBank accession numbers are 
shown in brackets. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 
replicates) is shown next to the branches. Phylogenetic analyses were conducted in MEGA4. The scale bar indicates nucle¬ 
otide substitutions per site. 
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FIG. 2. A comparison of part of the 3' untranslated region from three CoVs isolated from chickens (prefixed IBV), one CoV 
from turkeys (prefixed TCoV), and two CoVs from pheasants (prefixed ph), with the same region identified in the eight CoVs 
isolated from the captive birds described in this study. Dots indicate that all the viruses had the same nucleotide at that 
position. Deletions are marked with an asterisk. Yellow shading indicates the mutations observed in the Pavus/2/BR/ 
UNESP/09 sequence at a given position. 
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GU323510, Martinica/BR/UNESP/2009; GU323509, 
Martinica/BR/UNESP/2009. The CoV sequences AJ619580 
(pheasant CoV ph/UK/602/95), AJ619579 (pheasant CoV 
ph/UK/438/94), AJ310642 (Turkey CoV turkey/UK/412/ 
00), AJ278338 (IBV HVI-140), AJ278337 (IBV HV10), and 
AJ278335 (IBV D207) were obtained from GenBank for use in 
the phylogenetic analyses. Sequences were assembled using 
SeqMan Pro., and genome annotation was conducted with 
SeqBuilder (DNASTAR, Inc., v.8.0.2). Nucleotide and de¬ 
duced amino acid sequences were aligned using Clustal-W in 
the MegAlign program (DNASTAR, Inc.). A phylogenetic tree 
was constructed based on amino acid alignments using 
BLASTN analysis (www.blast.ncib.nlm.nih.gov/Blast.cgi). 

Results and Discussion 

Captive animals in zoos are part of a complex ecosystem 
that includes free-roaming animals and their associated in¬ 
fectious diseases, including respiratory viruses. These results, 
obtained from captive birds in Brazil, represent the first report 
of a circulating CoV inside a zoological park. The CoV RNA 
was detected in eight individual fecal samples (Table 1), 
which corresponds to an individual prevalence estimate of 
1.12% (95% confidence interval 0.5-2.7). Of the samples 
considered positive, two were from black shouldered pea¬ 
fowl (.Pavus cristatus), two were from American purple galli- 
nules (Porhyrula martinicd), and one each came from a 
Chestnut-bellied heron (Agamia agamia), an Emperor penguin 
(.Aptenodytes patagonicus), a Rusted margined guan (Penelope 
supercilliaris), and a Roadside hawk (Rupornis magnirostris). 
Interestingly, all of the positive samples were obtained from 
adult birds. Phylogenetic analyses were based on a sequence 
of 146 nt after the removal of primer sites (Fig. 1). Most re¬ 
markably, nine mutations were revealed in the Pavus/2/BR/ 
2009 (GU323507) sequence, despite many reports that con¬ 
sider the 3' UTR region to be highly conserved among group 3 
Co Vs (Cavanagh et al. 2002). One hypothesis that could ex¬ 
plain these mutations is the high density of different species of 
birds located in the same geographical area, which allows the 
virus to circulate among species (Jackwood et al. 2010). On the 
other hand, the natural habitat provides more distance be¬ 
tween species, particularly between predators and prey. In 
addition, the sequences of the CoVs isolated from the two 
peafowl samples clustered separately from the British de¬ 
posited sequences, as illustrated in Figure 1 (Cavanagh et al. 
2002). Sequences within this cluster were relatively homoge¬ 
nous between two European samples (Fig. 1), in accordance 
with previous studies (Cavanagh et al. 2002). 

Samples obtained from the Chestnut-bellied heron, two 
American purple gallinules, a Roadside hawk, and a rusted 
margined guan clustered together in a single clade, close to 
the sequence of the IBV D 207 strain (Figs. 1 and 2). In the 
present study of samples collected from 25 different species 
(kept in isolated locations within the zoological park), CoV 
RNA was detected in two Galliformes, two Ciconiiformes, 
one Gruiforme, and one Sphenisciforme. This demonstrates 
the high diversity in the susceptible species, as described 
previously (Woo et al. 2008). In fact, CoVs have been also 
detected in wildfowl species, rock doves, wild peafowl, 
waders, Amazon parrots, and some passerine species (Cava¬ 
nagh 2005, Liu et al. 2005, Gough et al. 2006, Woo et al. 2008, 


Hughes et al. 2009). The hypothesis that wild free-flying birds 
present in zoos may often act as vectors flying between en¬ 
closures that house the exotic birds and transmitting infec¬ 
tions is not fully understood. 

To understand the potential role of captive populations as 
reservoirs for CoVs, further studies are needed not only in 
South American zoo parks, but also worldwide. This is the 
first description of the molecular detection of group 3 CoVs 
related to captive birds in South America. The genetic diver¬ 
sity found in the black shouldered peafowl samples raises 
interesting questions about the 3' UTR conserved region, its 
evolution, and viral circulation among different species inside 
zoo parks. 
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